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Abstract 

We study the effect of a magnetic field on the behaviour of a conducting elastic rod subject to a 
novel set of boundary conditions that, in the case of a transversely isotropic rod, give rise to exact 
helical post-buckling solutions. The equations used are the geometrically exact Kirchhoff equations 
and both static (buckling) and dynamic (whirling) instability are considered. Critical loads are 
obtained explicitly and are given by a surprisingly simple formula. By solving the linearised 
equations about the (quasi-)stationary solutions we also find secondary instabilities described by 
(Hamiltonian-)Hopf bifurcations, the usual signature of incipient 'breathing' modes. The boundary 
conditions can also be used to generate and study helical solutions through traditional non-magnetic 
buckling due to compression, twist or whirl. 
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FIG. 1: Experimental setup for a conducting wire. 

I. INTRODUCTION 

A straight current-carrying wire held in tension between pole faces of a magnet is well 
known to buckle into a (roughly) helical configuration at a critical current (see Fig. [T]). A 
photograph of this phenomenon is shown in Section 10.4.3 of where a linear stability 
analysis is carried out for a simple string model. (A string is here meant to be a perfectly 
flexible elastic wire.) The problem was studied by Wolfe j^J by means of a rigorous bifurcation 
analysis for a (nonlinearly elastic) string suspended between fixed supports and placed in a 
uniform magnetic field directed parallel to the undeformed wire. He found that an infinite 
number of solution branches bifurcate from the trivial straight solution, much like in the 
Euler elastica under compressive load. In this case the non-trivial solutions are exact helices. 
That this should be so, is easily explained by the fact that the (Lorentz) body force is 
everywhere normal to the deformed configuration and hence the wire necessarily in a uniform 
state of tension. Some (statics) stability results (i.e., minimisation of the potential energy) 
were also obtained, indicating that the first branch of solutions is stable while the others are 
unstable. 



In a subsequent paper Wolfe {3] extends the analysis to a uniformly rotating (whirling) 
string and shows again the existence of bifurcating branches of whirling non-trivial solutions. 
Due to centrifugal effects no closed-form solutions could be obtained in this case. This result 
was further extended by Healey [4] using equivariant bifurcation theory in order to deal with 
the symmetries of the problem. 
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Wolfe also considered a conducting rod in a uniform magnetic field {5]. In addition to 
extension a rod can undergo flexure, torsion and shear, and for the case of welded boundary 
conditions it was found that in certain cases bifurcation occurs, with the usual infinity of 
non-trivial equilibrium states. All the works cited above were content with showing the 
existence of bifurcating solutions and did not study their post-buckling behaviour. 

In this paper we consider the post-buckling behaviour of a conducting rod. Wolfe consid- 
ered welded boundary conditions in J5] and encountered degeneracies (even-dimensional 
eigenspaces) because of rotational symmetry of the problem. In previous work [6] we 
showed that further complications occur and that magnetic buckling of a welded trans- 
versely isotropic rod (i.e., a rod with unequal bending stiffnesses about the two principal 
axes of its cross-section) is described by a remarkably degenerate pitchfork bifurcation. 
Wolfe also reported numerical evidence of helical post-buckling solutions. However, exact 
helical solutions cannot be supported by (coaxial) welded boundary conditions. Here we 
formulate a novel set of what we call 'coat hanger' boundary conditions that do support 
(i.e., are compatible with) exact helical solutions, and show that subject to these boundary 
conditions an isotropic rod does indeed buckle (exclusively) into a helix, or more precisely, 
that there is an infinite series of helical modes bifurcating at increasing load, each successive 
mode having one more (half) helical turn. All helical solutions can be obtained explicitly 
and it is found that the pitchfork bifurcations for these coat hanger boundary conditions are 
non-degenerate and that the critical loads are given by a remarkably simple formula. Unlike 
in string buckling a rod does not require a tensile force in the trivial state, but we allow for 
such an applied force as well. The pertinent dimensionless parameter that governs buckling 
measures the product of current and magnetic field against the bending force. 

We also study steady whirling solutions for which we introduce a rotating coordinate 
system. This extends Wolfe's analysis of whirling strings to whirling rods. An interesting 
feature of helical solutions is that since all points on a helix have equal distance to the 
whirling axis, and are therefore equally affected by centrifugal forces, solutions remain helical 
when spun. We perform a stability analysis by computing eigenvalues of the linearised 
boudary-value problem about a (quasi- stationary) whirling solution. For this we use a 
continuation (or homotopy) approach that takes advantage of the fact that exact expressions 
for the (imaginary) eigenvalues can be obtained in an appropriate limit (no spin, no magnetic 
field). The eigenvalues in this limit are then traced as system parameters are varied. 
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Whirl tends to destabilise the helical solutions, but stable solutions can be obtained by 
adding the effect of internal viscoelastic damping. We find Hopf bifurcations on the first 
bifurcating branch where a stable whirling solution becomes unstable under an increase 
of the angular velocity. We also briefly consider anisotropic rods. Critical loads can still 
be obtained analytically, but these rods buckle into coiled but non-helical solutions. Sec- 
ondary instabilities are found due to Hamiltonian-Hopf bifurcations, a common signature of 
'breathing' or 'flutter' instabilities in mechanical systems. 

Helical solutions are widely studied in a whole range of applications. Often these solutions 
are thought to arise through buckling of a straight rod under the action of end loads. 
However, as commented above, exact helical solutions are not supported by the usual set of 
boundary conditions. Consequently, boundary conditions are often not mentioned, or the 
rod is implicitly assumed to be infinitely long in order to prevent end effects [7] . One of the 
contributions of this paper is to present and highlight boundary conditions for an elastic rod 
that do support helical solutions. One could apply these experimentally if one was interested 
in generating or studying helical solutions in a finite-length rod. 

The paper is organised as follows. In Section 2 we present the rod mechanics formulation, 
in which the magnetic field enters the force balance equation through the Lorentz body force. 
The coat hanger boundary conditions are introduced together with a sketch of a device that 
can be constructed to realise these boundary conditions in a testing rig. For the study of 
whirling solutions the equilibrium equations are transformed to a coordinate system rotating 
at constant angular velocity. We use numerical bifurcation and continuation methods to find 
the buckling loads and to compute post-buckling solution paths, both for the statics and 
dynamics case. After the nondimensionalisation in Section 3, the linearisation is presented 
in Section 4. Section 5 first presents analytical stability results for the statics case and 
then introduces our continuation approach to numerical stability analysis of the full system. 
In Section 6 results are presented in the form of bifurcation diagrams and curves of Hopf 
bifurcations in appropriate parameter planes. Conclusions are drawn in Section 7 and the 
study closes with two Appendices: one giving details about the linearised system of equations 
and one deriving exact buckling results for helical solutions. The latter is complementary 
to the bifurcation analysis (not assuming any shape) in Section 5; together these analyses 
give a complete picture of helical magnetic buckling. 
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FIG. 2: Cosserat model of a rod. 



II. THE ROD MECHANICS MODEL 

We describe the elastic behaviour of a conducting cable by the Kirchhoff equations for 
the dynamics of thin rods. The rod is assumed to be uniform, inextensible, unshearable and 
intrinsically straight and prismatic. The assumptions of inextensibility and unshearability 
are appropriate for thin rods with relatively low external (here electroclynamic) forces. For 
the background of the Kirchhoff equations the reader is referred to 
were also used in a and 



These equations 



10[ to analyse the dynamics of a spinning tether. 
Let x denote the position of the rod's centreline and let {d 1? d 2 ,d 3 } be a right-handed 
orthonormal frame of directors (the Cosserat triad) defined at each point along the centreline. 
Since the centreline is assumed to be inextensible we can take ds in the direction of the local 
tangent: 

x'(s,t) = d 3 (s,t), (1) 

where the prime denotes differentiation with respect to arclength s measured along the 
centreline, and t is time. The directors di and d 2 will be taken to point along the principal 
bending axes of the cross-section (see Fig. [2]). The unstressed rod is taken to lie along the 
basis vector k of a fixed inertial frame {i, j, k}. 

Looking at Fig. [2] we note that the position vector of an arbitrary point of the rod can 
be expressed as 

X(s j ei,6 > *)=aj(s,t)+^idi(s,t)+6d 2 (s,t) 

(2) 

=x(s,t) + r(s, 

where (£1, £2) are the components of r in the cross-section relative to {di(s), d 2 (s)}. The rod 
is thus viewed as a set of infinitesimal slices centred at all s. A one-dimensional description 



will be obtained by averaging of forces and moments over each cross-section. The internal 
traction, which is the projection of the stress tensor onto the cross-sectional plane, is given 
by a force which we denote by / = /(s, £1, £2, (see Fig. [2]). The resultant elastic force 
exerted in a section S(s) is given by 

F(s,t)= [ f(s,^,^,t)dS, (3) 

JS(s) 

where dS is an infinitesimal area element. This force can be expressed in the director basis 
as F = Yli=i Q di. The resultant moment in the section S(s) is given by 

M(s, t)= [ r(s, £1, 6, t) x f(s, 6, £ 2 , t) dS, (4) 

JS(s) 

and will be expressed as M = Y2^=i Midi. 

The rod is assumed to carry an electric current for which we can write 

I = Ix' = Id 3 . (5) 

Here we have assumed the current to have the same direction as the rod, which is consistent 
with a one-dimensional rod theory. It amounts to the assumption that the cross-section 
of the conducting wire is small enough to make currents within the cross-section (eddy 
currents) induced by the motion negligible. The current J interacts with the magnetic field 
-Bo to generate a (Lorentz) body force given by 

F L = Id 3 x B (6) 

Following Q we assume the magnetic field to be uniform and directed along the unstressed 
rod, i.e., 

B = B k. (7) 

The balancing of forces and moments across an infinitesimal rod element then yields the 
following set of partial differential equations js, 9|: 

F 1 + I Bods xk = pAx, (8) 

M' + d 3 x F = p{I 2 d x x d\ + hd 2 x dk), (9) 

where p is the (volumetric) mass density, A the cross-sectional area, Ii and 7 2 the second mo- 
ment of area of the cross-section about d\ and d 2 respectively, and ( ) denotes differentiation 
with respect to time. 
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For a closed system of equations these balance equations need to be supplemented by 
constitutive relations that characterise the material behaviour of the rod. We assume the 
rod to be made of homogeneous isotropic linear viscoelastic material so that stress-strain 
relations, based on a model by Valverde et al. [10], are 

Mt = £Ji(«i+7 w /ci), 

M 2 = EI 2 (k 2 + 1v k 2 ), (10) 
M 3 = GJ(k 3 + 7„k 3 ), 

where K\ and k 2 are the curvatures about d\ and d 2 , respectively, while k 3 is the twist about 
g?3. The constant 7^ is the viscoelastic coefficient of the material, E is Young's modulus, 
G is the shear modulus and J is the second moment of area of the section about d 3 . We 
shall assume that the section is symmetric with respect to the principal axes, in which case 
J = h + h- 

The Ki £1X6 the components of the curvature vector 

3 

K = y^Kj dj, (11) 

which governs the evolution in space of the frame of directors as one moves along the 
centreline: 

d'^KXdi (i = 1,2,3). (12) 

The constitutive relations ffTO]) can be used to replace the k« in (T121 by moments, after 
which the equations (pQ), (jHJ), (jSJ) and (|T2|) form a system of 18 differential equations for the 
18 unknowns (x, F, M, d\,d 2 , d 3 ). 



Remark: We ignore in this study secondary electrodynamic effects (such as an induced emf 
and hence additional current in the conductor) as a result of the motion of the wire in the 
magnetic field [TJ|. Since we are considering a steadily rotating wire these effects would be 
null on the configuration of the wire. However, the same would not be true for the stability 
analysis, which considers arbitrary time- dependent perturbations. We assume that these 
induction effects are negligible. 
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A. Equations of motion in a uniformly rotating frame 



We shall also be interested in steadily rotating solutions and therefore we transform 
the equilibrium equations (jSJ) and OH]) to a coordinate frame {e 1; e 2 ,e 3 } that rotates with 
constant angular velocity uj = uk about the k axis (and the axis of the rod in its trivial 
unstressed state). Noting that the derivative with respect to time of an arbitrary vector 
V(s, t) is given by 

+ uxV(s,t), (13) 



dV{s,t) 



dt 



dV(s,t) 



dt 

where 4 1 . indicates the derivative with respect to time in the inertial frame and 4 1 stands 

at \ i 1 dt\m 

for the derivative with respect to time in the moving frame, the equations (|S]) and (JH]) 
expressed relative to {ei,e 2 ,e 3 } become 

F' + IB d 3 x e 3 = pA(x + 2uj x x + oj x (u> x x)), (14) 
M' + d 3 x F = ph{d\ x d\ + 2di x (w x dt) + (w • x u>)) 

(15) 

+ p/i(d 2 x d 2 + 2d 2 x (w x da) + (w • d 2 )(d2 x w)). 
The second term on the right-hand side of equation (fl4"l) is the Coriolis force, while the third 
term is the centrifugal force, as a result of the rotating coordinate system. 

Steadily rotating (whirling) solutions satisfy the equations ( IT41 and (TT5*1) with the dotted 
variables set to zero: 

F' + IB d 3 x e 3 = pAw x (a; x a;), (16) 

M' + d 3 x F = p/ 2 (w • dt)(dt x u>) + p/i(a> ■ d 2 )(d 2 x a>). (17) 

The other equations jT]) and ( Fl2i) do not change their form, but all vectors are now to be 
considered as expressed relative to the rotating frame {ei,e 2 ,e 3 }. Statical solutions are 
simply obtained by setting u> equal to zero. 

For a well-posed problem the final 18 ODEs require 18 boundary conditions to be speci- 
fied, which we do next. 

B. Coat hanger boundary conditions 

Helical solutions in rods are usually studied in infinitely long rods, which avoids the need 
for imposing boundary conditions. Indeed, it is not immediately clear how an exact helix 
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FIG. 3: Coat hanger boundary conditions. 

can be supported: the boundary conditions cannot be simply welded as no two points on 
a helix have coaxial tangents, nor can they be simply pinned because a helix has curvature 
and therefore carries a bending moment. Here we formulate a set of boundary conditions 
that support exact helical solutions. We call them coat hanger boundary conditions, for 
obvious reasons. 

Consider Fig. [3] where a rod is suspended between two axes v Q (at s = 0) and V\ (at 
s — 1) lying in two parallel planes normal to e 3 . Axis v% is taken to be fixed in space, 
while Vq is free to move along e$. We assume the axes to have a fixed relative rotation x, 
i.e., Vq ■ V\ = cosx- The rod is free to hinge about and slide along both Vq and v±. For 
definiteness we assume that the rod is mounted in such a way that both axes Vq and v\ are 
directed along the vector d 2 in the rod's cross-section. This situation is described by the 
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following boundary conditions: 

di(0,f).«o = 0, (18) 

d 3 (0,t)-v = 0, (19) 

M(0,t)-v = 0, (20) 

F(0,t)-v = 0, (21) 

F(0,t)-e 3 = T, (22) 

x(0,t)-(« o x e 3 ) = 0, (23) 

at s — 0, and 

di(L,*)-«i = 0, (24) 

d3(L,t)-«i = 0, (25) 

M(L,t).«i = 0, (26) 

F(L,t). Vl = 0, (27) 

x(L,t) ■ (vi x e 3 ) = 0, (28) 

z(L,t) = L, (29) 

at s = 1, where the position vector has been decomposed as x = xe.\ + ye 2 + ze% and T is an 
applied end force (positive for tension). Conditions (1231) and (|28|) restrict the movement of 
the ends of the rod to the planes spanned by (v , e 3 ) and (i>i, e 3 ), respectively. To these 12 
conditions we have to add conditions that ensure the orthonormality of the director basis, 
for which we can take 

di(0,*)-di(0,*) = 1, 
d 2 (0,t) ■ d 2 (0,t) = 1, 
d 3 (0,t)-d 3 (0,t) = l, 

(30) 

di(0,t) ■ d 2 (0,t) = 0, 

di(0,t)-d 3 (0,t) = 0, 

d 2 (0,t)-d 3 (0,t) = 0, 

for a total of 18 boundary conditions, as required. 

We shall take x — so that the initial rod, lying straight along e 3 , is untwisted. This 
choice implies that any helical solutions will have an integer number of half helical periods. 
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Also, the directors, and hence the cross-section of the rod, will make a half-integer number 
of turns between s = and s = L. Without loss of generality we may choose Vo — V\ — e 2 , 
so that at the ends of the initial rod the directors {di, d 2 , d 3 } are aligned with {ei, e 2 , 63}. 

Since v and Vi are chosen aligned, there is a rigid-body degree of freedom of translation 
of any solution along these axes. To eliminate this degeneracy we replace condition (J27J) 
above by 

x(L,t) -Vi = n, (31) 
where the slide r\ along V\ (see Fig. [3]) is chosen as follows. First note that a helix is a curve 



of constant axial radius r = ^ x 2 + y 2 and (total) curvature k = \f k\ + k\. The two are 
related by = sin 2 9, where 9 is the helical angle defined by d^-e^ = cos 9 if is along the 
axis of the helix (the angle | — 9 is usually called the pitch angle). Since the rod is hinged 
about d 2 , we have n 2 (L,t) = 0. So, in order to ensure that any bifurcating helix is centred 
at e 3 we take 

n = '-WM->' , (32) 

Ki{L,t) 

giving a nonlinear boundary condition. Note that K\ can here be taken with its sign, so that 
( 13T1) also specifies which way the rod moves along V\. When the rod buckles, the axis 
lifts up and the rod is free to find its own radius r = |ri|. We stress that condition (|3TT) has 
no effect on the bifurcation behaviour. In particular, it does not suppress any non-helical 
solutions. It merely ensures that if a helical solution bifurcates it will be centred at the 
axis of rotation. This is important when we start rotating the axes v$ and v\ about e%. A 
centred helix will experience a uniform centrifugal force and is therefore expected to remain 
helical. 

Of course the above coat hanger boundary conditions merely allow for helical solutions. 
They need not exist. However, if the equilibrium equations do have helical solutions and 
one-parameter curves of such solutions intersect the trivial path of straight solutions, then 
one might expect to detect them as (pitchfork) bifurcations at critical buckling loads. The 
results presented in Section |VT] show that this is indeed the case. 
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III. NONDIMENSIONALISATION 



We make the system of equations dimensionless by scaling the variables in the following 

way 



, Eli - s r , _ u _ x 



(33) 



Here u c is a reference characteristic bending frequency of the rod. 

With this nondimensionalisation the equations become (dropping the overbars for sim- 
plicity and letting a prime denote j| and an overdot 



F' + Bd 3 x k = x + 2lj x x + w x (u> x x), (34) 



M' + d 3 x F = P R(di x d\ + 2d 1 x u> x <ii + (a? • di)(di x a;)) 

+ (d 2 x d 2 + 2d 2 x a? x d 2 + (w • d2)(d 2 x u>)) , (35) 

= da, (36) 

d' i = kx d u (37) 

and the constitutive relations can be written as 

r(i + R) 

M = («! + 7 «i)di + i?(« 2 + 7^ 2 )d 2 + H« 3 + 7^s)d 3 , (38) 

where the dimensionless parameters are 

h h n B IL* 2G 

P= AL^ R= V B = ^h' r = ^' 7 = W ' (39) 

and — 1) is equal to Poisson's ratio. For the boundary conditions we can still use ([18]) 
to (1291) if we assume that they now refer to dimensionless variables and that the right-hand 
conditions are imposed at s = 1. 



IV. PERTURBATION SCHEME LINEARISATION 

We consider whirling solutions (relative equilibria) that are stationary in the moving 
frame {ei, e 2 , e 3 }. Such solutions are found by solving the set of equations (I34"l) - (l38"]) with 
the dotted variables set to zero (thus obtaining an ODE). To study their stability we linearise 
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hese whirling solutions. We follow the approach in 10], which 



the full PDE dMD-(|SHD about 
is similar to approaches in 7, [12|. The stability of static (non- whirling) solutions can be 
investigated by simply setting the angular velocity uj to zero. 
We start our perturbation analysis by writing 



di(s,t) = d°{s) + 5d t l {s,t) + 0(5 2 



i = 1,2,3, 



(40) 



where d®(s) are the basis vectors of a quasi-stationary solution, d\(s, t) are the basis vectors 
of a time-dependent perturbation and 5 is a small bookkeeping parameter introduced to 
separate scales. Note that, in order to preserve orthonormality to 0(5) (di-dj = 5^ + 0(5 2 )), 
we must have 



dl(s,t) = J2AMt)d%s) 



* = 1,2,3, 



(41) 



where the matrix Ay is skew-symmetric and can be written as 

^ a 3 — a 2 ^ 



V 



(42) 



— «3 «i 
«2 —OL\ 

Thus, the nine components of the director basis perturbation are described by only three 
independent parameters, and if we introduce 



a 



( a u a 2 , a 3 ) 



(43) 



(with respect to the unperturbed director basis) then the perturbed director basis can be 
expressed as 

di(s,t) = d° i {s) + 5a{s,t) x d°(s) + 0(<5 2 ), i = 1,2,3. (44) 

Using (|44|) . the perturbation of an arbitrary vector V = Y^=i Vidi can be written on the 
basis {d^d^cZ^} as 



V = V° + 5V l + 0(5 2 ) = J2lK° + SO*? + (« x V\M + 0(5 2 ), 



(45) 



i=l 



where ()j denotes the component along d® and time and space dependence of the variables 
have been suppressed for the sake of simplicity [lo1 |. 

Applying this perturbation scheme to the PDEs (13"4"l) - (f3"5l) and the boundary conditions, 
we arrive at an 0(1) nonlinear ODE for the quasi-stationary solutions and an 0(5) linear 
PDE governing their stability. 
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A. The O(l) equations — quasi-stationary whirl 



The 0(1) equations are time- independent. Recalling that u) = ue^, we find the 0(1) 
terms of the linear momentum equation ( 13~4"|) . projected on the director basis {d®, d^, d 3 }, 
to give 

{FfY - + F°4 + B(d° 32 d° u - d 31 d° 12 ) = ~co 2 (x d° u + y°d° 12 ), (46) 
(F 2 )' - F 3 ° K ? + F?k° 3 + B(d° 32 d° 21 - d° 31 d° 22 ) = -co 2 (x°d 21 + y°d° 2 2), (47) 
(F 3 )' - F?4 + F$4 = -co 2 (x°d° u + yX), (48) 

where subscripts are used to indicate components relative to the basis vectors {d 1 ,d 2 ,d 3 y } 
(but the d® components are relative to {e 1; e 2 , e 3 }). Similarly, the 0(1) term of the angular 
momentum equation (13"51) . projected on the director basis {d®, d 2 , d 3 } gives 



2 M M M° M° 

M)' = j^j- - ^ + F 2 ° + P^dlMA, - dU° 12 ), (49) 

{M ° y = -f^TR) + M ° M ° - F " + Pi ^ 2 <(« - (so) 

{Ml)' = - M°M$ + PRoo^MA, - d° u d° 3 2) + P^dlMA, - d° 21 d°32)- (51) 



The 0(1) term of equation ( 1361) can be expressed as 

(*°)' = 4 (52) 

and the twist equation (13"T|) by 

(d°)' = K °xd°, i = 1,2,3, (53) 
where k, = Y^j=i rfd®. The 0(1) term of constitutive relations ( 138]) can be expressed as 

M° = < M2 = Rk», M° 3 =^±^4, (54) 
which can be used to express the k° in (1531) in terms of the moments M 4 °. 
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Proceeding in the same way, the 0(1) part of the boundary conditions is given by 



4(0) ■ v = 0, (55) 

4(0) ■ v = 0, (56) 

M°(0) • v = 0, (57) 

F°(0) ■ v = 0, (58) 

F°(0) • e 3 = T, (59) 

x°(0) ■ (v x e 3 ) = 0, (60) 

= 0, (61) 

d°(l) • v x = 0, (62) 

M°(l) • «! = 0, (63) 

x°(l) • ( Vl x e 3 ) = 0, (64) 

x°(l)- Vl = r°(l) = -i-(l-(rf° 3 (l)) 2 ), (65) 

z°(l) = 1, (66) 



where (|65ll is the 0(1) contribution from ( 1311) . with the radius of the helix also affected by 
the perturbation scheme, i.e., r(s,t) = r°(s) + <5r'(s,t). 

B. The 0(5) equations — linearisation 

The 0(5) part of the linear momentum equation (|34p can be written as 

(F\s, t))'+B 1 {8)F t (8, t) + B 2 (s)x t (s, t) + B 3 (s)a'(s, t) + B 4 (s)cx(s, t) 

(67) 

=B 5 (s)x t (s,t) + B 6 (s)i*(s,t), 

where the 3x3 matrices B^(s) are given in Appendix A. Here we have expressed F l relative 
to {4-, 4i 4} an d xt relative to {ei, e 3 }. For the 0(8) part of the angular momentum 
equation (|35l) we can write 

(M t (s,t))'+C 1 (s)M t (s,t) + C 2 (s)a'(s,t) + C 3 (s)a(s,t) + C 4 (s)F*(s,t) 

(68) 

=C 5 (s)a(s,t) + C 6 (s)a(s,t), 
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where the matrices Cj(s) are again given in Appendix A. M t is expressed relative to 
{d\, e?2, e^}. The 9 twist equations (l37j) at 0(6) are reduced to only 3 independent equations 
that relate re* and a as 

K*(a, t) = a'(s, *) + K °( s ) x a ( s i t). (69) 
Introducing these relations into the 0(6) part of the constitutive relations gives 

M\s, t) + D 1 (s)a'(s 1 1) + D 2 (s)a(s, t) = D 3 (s)a(s, t) + D 4 (s) — (a(s, t)), (70) 

where the matrices Di(s) are given in Appendix A. Finally, the 0(6) part of equation ( 1361) 
yields 



(**(*,*))' = a(M)xd°( a ). (71) 
Applying the perturbation scheme to the boundary conditions at 0(6), we obtain 



y\^t) + 



(dg(l).cs)(4(l,*).cs) + 



ai(0, t) 


= 0, 


(72) 


a 3 (0,i) 


= 0, 


(73) 


x*(0,t) 


= o, 


(74) 


)a 2 (0,t) 


= 0, 


(75) 


F 2 '(0,t) 


= 0, 


(76) 


M*(0,t) 


= 0, 


(77) 


z*(l,t) 


= 0, 


(78) 


z\l,t) 


= 0, 


(79) 


ai(l, i) 


= 0, 


(80) 


a 3 (l,t) 


= 0, 


(81) 


M*(l,t) 


= 0, 


(82) 


)«*(!,*) 


= 0. 


(83) 



After elimination of the k° by means of (l54l . the set of 12 equations (IBTl) . ( J68l) . ( 1701) and 
( TTTT) together with the 12 boundary conditions (1721) -(1531. with appropriate initial conditions 
form a well-posed initial-boundary-value problem. 
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V. STABILITY ANALYSIS 



Since we are interested in stability of solutions we look for solutions of the 0(5) equations 
of the form 



When these expressions are inserted into (j6"Tj) - (l7Tl) a linear eigenvalue problem for a 12- 
dimensional ODE is obtained in terms of the variables (x l , a*, F ,M ). The eigenvalue 
A measures the growth of small perturbations and is to be found as part of the solution. 
Eigenvalues come as complex conjugate pairs. A whirling state is unstable if at least one of 
the (in general infinitely many) A's has positive real part. 

To solve a real system of equations we split the eigenvalues and variables (eigenfunctions) 
into real and imaginary parts, A = A r + x l = x l r + ix\, a 1 = a* + idt\, F = F r + %F i 
and M = M r + iM { . The equations (1671)-(I7T|) along with the boundary conditions (1721)- 
( |83l) . are similarly split into real and imaginary parts. Thus we end up with a doubled 
24-dimensional linearised boundary-value problem. 

A. Stability of the straight rod — static magnetic buckling 

The trivial solution of the 0(1) equations ( |46|) -( |5^|) . representing a straight and untwisted 
rod, is given by 

x{s) = se 3 , F(s) = -Te 3 , M(s) = 0, d^s) = a (i = 1,2,3), se[0,l). (88) 

It satisfies the coat hanger boundary conditions. For the statics case (u = 0) without end 
force (T = 0) the 0(5) equations f l6T|) . f lBS]) . f lTU|) . fl7T|) about this trivial solution, on inserting 
(18^1) -(IHTj). can be written as 



x\s,t) 



= x (s)e 



xt 



(84) 



et(s,t) = et(s)e , 
F t (s,t) = F\s)e xt , 
M\s,t) = M(s)e xt . 



(85) 
(86) 
(87) 




Ml - —M 3 = 0, 



(89) 
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with boundary conditions 

x (0) = x(l) = x"(0) = x"(l) = y"'(0) = y(l) = y'(0) = y'(l) = 0, M^O) = = AfJ(l), 

(90) 

while F 3 = 0, z = 0. Note that the torsional (M 3 ) modes decouple from the bending 
modes. 

To find the static magnetic buckling loads we set A = 0. The bending equations then 
reduce to 

z""" + ^-z = 0, for z = x', (91) 
R 

subject to 

x (0) = x(l) = x"(0) = x"(l) = x""(0) = x""{\) = 0. (92) 

On setting z = e tks we obtain the characteristic equation — k 6 + B 2 /R = with solutions 
k lj2 = ±p, ^3,4,5,6 = ±P (l ±i>/3) /2, where /3 = B 1 / 3 /R 1 / 6 . Application of the boundary 
conditions (192]) to the general solution z(s) = Y^j=i a j e%kjS leads to the remarkably simple 
condition: 

B = (nn) 3 VR. (93) 

These critical loads correspond to pitchfork bifurcations where non-trivial solutions bifurcate 
from the trivial straight solution. We stress that the above calculation is only possible for 
the statics case. If u ^ then the x and y equations do not decouple and no simple 
characteristic equation is obtained. However, bifurcating branches of helical solutions, and 
hence critical loads, can be computed explicitly, even for non-zero u; see Appendix B. 

B. Eigenvalues for the unperturbed problem (T = 0, 7 = 0, uj = 0, B = 0) 

We shall call the case where T = 0, 7 = 0, u = and B = the unperturbed problem. 
For this problem explicit expressions can be obtained for the eigenvalues of the linearisation 
about the straight solution. The x and y equations in (|89l) decouple into two fourth-order 
beam equations: 

x"" - X 2 Px" + ^-x = 0, 

R (94) 

y"" - \ 2 Py" + \ 2 y = 0, 
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subject to boundary conditions ( 190]) . Since we anticipate imaginary eigenvalues we set 
A = ifi, x = e tks , y = e tKS , and find for the x equation 

1/2 



/l 1 \ 1 

k h2 = ± i-fx 2 P + -v//i 4 P 2 + V/i?] =: ±a, 

fX 1 \ 1/2 

k 3A = ±i l-^^P 2 + 4^/R--^P\ =: ±ib 



while for the y equation 

1/2 



Kl52 = ± i-fi 2 P + -v^AP^+VJ =: ±«, 

/I X \l/2 

where a, b, a, (3 are non-negative real numbers. The general solutions are 

x(s) = A x sin as + B x cos as + sinh 6s + D x cosh 6s, 
y(s) = A y sin as + B y cos as + C y sinh (3s + D y cosh (3s. 

Application of the boundary conditions fl90|) leads to 

sin a = and (3 cos a sinh /3 + a sin a cosh/3 = 0. 

The first equation implies 



n 2 7r 2 



/^ = ±^^= 9 „ 1/r> > n = 1,2,3, .. 
'nH 2 P + 1/R 



The second equation is transcendental and needs to be solved numerically to obtain the 
eigenvalues /i, for instance by using a Newton-Raphson scheme. Meanwhile, the torsional 
eigenvalues for the M3 equation in ( 1891) are given by 



, r 

H = ±nix\(— , n=l,2,3, 



These are all the eigenvalues for the unperturbed problem. They will be used as starting 
values in the numerical procedure described next. 



C. Numerical procedure 



The main idea is to use the known eigenvalues in the unperturbed problem as starting 
values in a continuation procedure in order to compute the eigenvalues and corresponding 
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eigenf unctions for general values of the parameters T, 7, 00 and B. For this we use the well- 
tested code AUTO [13| (specifically AUTO2000). AUTO solves boundary- value problems 
by means of orthogonal collocation. It requires a starting solution and can then trace out 
solution curves as a parameter of the problem is varied. Bifurcations are detected where 
branches of solutions intersect. At such points AUTO is able to switch branches and compute 
curves of bifurcating solutions. 

Our procedure takes advantage of the fact that A appears only quadratically in the 
linearisation (167|) . (|68|) . (170|) and (ITTj) if 7 = and u = 0. To explain the method consider 
the typical 0(5) equation 

z"" - X 2 f(s)z" + X 2 g(s)z = 0, (95) 

where / and g are functions of the 0(1) solution. Writing z = x + iy, A = A r + zAj, we can 
decompose the z equation into 

x"" - (A r 2 - X 2 )f(s)x" + 2X r X l f{s)y" + (A r 2 - X 2 )g(s)x - 2X r X ig (s)y = 0, 

(96) 

y »» _ ( A 2 _ \*)f( s ) y " - 2X r X t f(s)x" + (A r 2 - X 2 )g(s)y + 2X r X ig (s)x = 0. 
The important thing to note here is that these equations decouple into two identical equa- 
tions if the eigenvalue is either imaginary (A. r = 0) or real (Aj = 0). 

This suggests the following sequence of steps, involving boundary-value problem of in- 
creasing dimension, to compute eigenvalues of static or uniformly whirling solutions. 

1. Consider the unperturbed problem of Section TV Bl and, noting that all eigenvalues are 
purely imaginary, solve the 30-dimensional system of 18 0(1) equations and one 12- 
dimensional system for the imaginary part of the 0(5) equations (cf. the y equation in 
(|96|) ). Set A r = and use Aj as the continuation parameter in AUTO in order to com- 
pute the eigenvalues (instead of solving the transcendental equations in Section l"VBj) . 
These eigenvalues will show up as branching points (BP), or pitchfork bifurcations, as 
eigenvalues by definition are those values for which non-zero BVP solutions exist. By 
symmetry it is only necessary to consider Aj > 0. 

2. Keeping the same 30-dimensional system, switch branches at a BP to compute ('grow') 
the corresponding (imaginary) eigenfunction. Since the equations are linear the value 
of Aj will not change in this run. For later use we monitor the non-zero solu- 
tion by means of some measure ||.||j (not necesarily a proper norm) on the space 
{a;-, 6t\, .Fj, M { } of imaginary linearised variables. 
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3. Now consider the full system of 42 equations (18 0(1) equations and two sets of 12- 
dimensional 0(5) equations (cf. (196]) ). Fix the measure ||.||j on the imaginary part 
and release X r instead in order to compute the real eigenfunction (since the imaginary 
part of the solution is fixed there is only one branch of solutions through the starting 
point and there is nowhere else to go for the continuation but to 'grow' the real 
eigenfunction). Again we monitor this function by means of a suitable measure ||.|| r - 
In this run neither A r nor Aj will change. 

This approach works because the solution obtained in step 2 also solves the full 42- 
dimensional system when the extra 12 variables (x l r , ck*, F r , M r ) are set to zero. This 
is a consequence of the fact that the real and imaginary parts of the 0(5) equations 
decouple if A r = 0, as a result of the quadratic dependence of the eigenvalue problem 
on A (cf. (JHSD). 

Steps 2 and 3 can be performed for as many of the BPs computed in step 1 as required and 
will give the corresponding eigenvalues and eigenfunctions. Once these have been obtained 
both measures ||.||. r and ||.||j can be fixed and an extra system parameter such as B or u 
released in order to trace the eigenvalues (and hence monitor stability changes) as system 
parameters are varied. (Note that fixing ||.|| r and ||.||j makes sense as eigenfunctions are 
only defined up to a multiplicative factor.) 

The above 3-step procedure is not limited to linearisations about the trivial straight 
solution. It can be applied to any starting solution that has no eigenvalue with both A r 
and Aj non-zero, as these would not be picked up in step 1. (It is of course no problem if 
eigenvalues become fully complex (Hopf bifurcation) in the course of further continuations.) 
For instance, we find that at the first critical B, given by (l9~3"|) . the lowest conjugate pair of 
eigenvalues ±Aj goes to zero and becomes a real pair of eigenvalues, signalling a stability 
change of the straight rod. The (first-mode) solution bifurcating at this point is stable 
with all eigenvalues being imaginary and the above procedure can be applied to find the 
eigenvalues. 

We end this section with a few comments: 

(i) There are infinitely many eigenvalues and the above procedure only finds the lowest 
order ones. This is of course a limitation of any numerical scheme. We find that 
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FIG. 4: Bifurcation diagram for a stationary isotropic rod subject to coat hanger boundary condi- 
tions (cj = 0). 

eigenvalues vary slowly with system parameters, suggesting that stability is governed 
by the lowest-order eigenvalues. We typically consider 5 or 6 eigenvalues. 

(ii) Note that in steps 1 and 2 above we could not have taken the full 42-dimensional system 
of equations as that would have made the branching points (pitchfork bifurcations) 
degenerate and AUTO would not detect a BP. This is because if A r = (or Aj = 0) 
the two sets of 12-dimensional linearised equations are identical (cf. (196]) ). 

VI. NUMERICAL RESULTS 

A. The statics case (oj = 0) 

Fig. H] shows the bifurcation diagram obtained by varying the magnetic field parameter B 
in the isotropic case (R = 1). Throughout this entire section the dimensionless parameters 
taken, unless stated otherwise, are those listed in Table HI where realistic dimensional pa- 
rameters are also given. Pitchfork bifurcations on the trivial branch occur at BP1 = 31.01, 
BP2 = 248.05, BP3 = 837.17, etc., in agreement with (193]) (only the positive-r branch is 
shown). The bifurcating solutions are found to be exact helices and therefore the dimension- 
less helical radius r/ Lis used as solution measure on the vertical axis. A slight complication 
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TABLE I: Parameters used for the coat hanger boundary conditions. 



L 


5 m 


P 


0.001 


A 


3 x 1(T 5 m 2 


R 


1 


E 


30 x 10 9 N/m 2 


r 


0.76923 


h = h 


9 x KT 11 m 4 








x y 



FIG. 5: The {ei-e3} and {e2-e3} projections of the first 4 helical modes at constant curvature 
k = 8. Values of B are: 61.68 (bl), 270.69 (b2), 868.72 (b3) and 2024.70 (b4). 

in computing this diagram occurs because of the denominator in ( 1321) . which is zero for 
the straight rod. However, this problem is easily resolved by replacing boundary condition 
( |311 by x(L) ■ vi = along the trivial branch and switching back to ( |311 once an incipient 
non-trivial solution has been obtained. 

Fig. [5] shows {ei-e^} and {e2-e%} projections of bifurcating solutions along the first four 
branches, taken at constant curvature k = 8. It was noted in Section III Bl that the coat 
hanger boundary conditions with x — (i.e., parallel end supports v and only allow 
helices of a half-integer number, n, of helical turns. We find that each successive bifurcating 
solution in Fig. |4] has one more half helical turn. For large B the solutions approach a 
circular shape in the z = 1 plane, with corresponding value r/L = l/(n7r) along the vertical 
axis in Fig. HI The bifurcating branches have handedness. That is, regardless of the sign 
of r, the bifurcating solutions are right-handed helices if, as here, B > and would be 
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FIG. 7: Evolution of imaginary (a) and real (b) parts of the first five pairs of eigenvalues along the 
trivial solution, from B = to B > BP2. 

left-handed helices if B < 0, i.e., had we run B in the other direction. Fig. [6] shows three- 
dimensional views of two solutions along the sixth branch. They have three full turns and 
nicely illustrate the exact helical shape. 

Fig. [7] show the evolution of the imaginary and real parts, respectively, of the first five 
pairs of eigenvalues along the trivial solution, from B = to B > BP2. At successive 
pitchfork bifurcations, pairs of imaginary eigenvalues collide at zero and become real, one of 
the eigenvalues of the pair with positive real part, signalling that the trivial solution becomes 
unstable at B = BP1, while a further loss of stability occurs at BP2 where a second pair 
of eigenvalues becomes real. 

Fig. [H] shows the evolution of the first five pairs of imaginary eigenvalues when switching 
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FIG. 8: Imaginary part evolution of the first five pairs of eigenvalues along the trivial solution, 
from B = to B = BP1 and then switching to branch bl. All real parts of the eigenvalues are 
zero. 
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FIG. 9: Evolution of the imaginary (a) and (positive) real (b) parts of the first five pairs of 
eigenvalues along branch b2. 

at BP1 from the trivial branch to bl, which is found to be stable. As expected, solutions 
along b2 are found to be unstable (see Figs[9k,b). 

B. Whirling solutions (w 7^ 0) - Hopf bifurcations 

Next we consider whirling motions in which the supports Vq and V\ are spun about k 
with constant angular velocity u. Solutions remain helical, as expected. Fig. HO] shows the 
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FIG. 10: Effect of u on the helical radius of the solution along bl in Fig. [5] at B = 35. The dot 
indicates a Hopf bifurcation where the real parts of two complex conjugate eigenvalues become 
positive (see Fig. PHT a)) and the solution loses stability. 

effect of to on a solution taken on the first bifurcating branch (bl) of Fig. HI The helical 
radius increases with u. Fig. [11] shows the bifurcation diagram for a fixed value of u = 2. 
Pitchfork bifurcations along the trivial branch occur at B = 29.75, 247.44, 836.80, etc. 

The angular velocity tends to destabilise the helical solutions (see Fig. [12], where the 
eigenvalues are shown as a function of w). We now introduce damping by continuing in the 
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FIG. 11: Bifurcation diagram for a whirling isotropic rod subject to coat hanger boundary condi- 
tions (oj = 2). 
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FIG. 12: Evolution of imaginary (a) and real part (b-f) of the first five pairs of eigenvalues with 
respect to uj for a solution of the first branch bl of Fig. 2] for B = 35, bifurcation diagram in 
Fig. QUI Note that first panel only represents the positive imaginary eigenvalue, as the negative is 
symmetric with respect to the real axis. 
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FIG. 13: Evolution with respect to u of the imaginary part of the first five pairs of eigenvalues for 
a solution along bl at B = 35 and 7 = 0.05. The first two panels only show the positive imaginary 
part. 

parameter 7, and to further investigate the stability of the first mode we fix 7 = 0.05 [h]] 
and perform continuation in u. Figs [13] and [TH show that at u = the solution is stable 
and that around u = 1.75 the real part of the first eigenvalue becomes positive. Thus the 
system loses stability in a Hopf bifurcation. The point where this occurs is indicated in the 
bifurcation diagram in Fig. [TUJ which is still valid as 7 has no effect on relative equilibria. 
Fig. [15] shows curves of Hopf bifurcations in the co-B parameter plane for various values of 
7. At the end points of these curves (indicated by dots) the Hopf bifurcation coalesces with 
the pitchfork bifurcation in which the curve bl is created. Consequently, the whirling helical 
solution is stable below these curves and unstable above. For small 7 the curve collapses 
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FIG. 14: Evolution with respect to u of the real part of the first five pairs of eigenvalues for a 
solution along bl at B = 35 and 7 = 0.05. 

onto the line u = 0. Note that for relatively large u damping has a stabilising effect, but 



that for small positive and for negative u damping has a destabilising effect. The 



behaviour is known from the classical linear stability theory of gyroscopic systems 



after 



C. The stationary anisotropic rod — Hamiltonian-Hopf bifurcations 

It is known that in the absence of damping and inertial effects (i.e., 7 = 0, oj — 0) the 

15j. In Hamiltonian 



equations for a rod in a magnetic field have a Hamiltonian structure 
systems a common mechanism for loss of stability is through a so-called Hamiltonian-Hopf 
bifurcation 16j. In this bifurcation two imaginary eigenvalues move on the imaginary axis, 



meet at some non-zero value and then leave the axis to become a real pair, as illustrated 
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FIG. 15: Loci of Hopf bifurcations along bl for different values of the damping parameter 7: from 
left to right, 7 = 0, 7 = 0.01, 7 = 0.025 and 7 = 0.05. Dots indicate points where the Hopf 
bifurcation coalesces with the primary pitchfork bifurcation. 
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FIG. 16: Eigenvalue behaviour in a Hamiltonian-Hopf bifurcation. 

in Fig. [16] Since eigenvalues come as conjugate pairs this event involves four eigenvalues. 
In structural problems Hamiltonian-Hopf bifurcations usually mark oscillatory instabilities 
such as flutter. 

No Hamiltonian-Hopf bifurcation was found in our study of the isotropic rod in Section 
IVI Al However, after breaking the cross-sectional symmetry by taking R < 1 Hamiltonian- 
Hopf bifurcations are found to occur under B continuation. Fig. [17] shows a bifurcation 
diagram for R = 0.5512. Pitchfork bifurcations along the trivial branch are found at B = 
23.058, 184.56, 622.55, etc., in agreement with (|93|) . Bifurcating solutions are not pure 
helices. Figs [TH] and [TH] show the eigenvalue behaviour for a solution on the first bifurcating 
branch, bl, in Fig. [T7J The first and second eigenvalues collide at B — 23.26, where the 
solution becomes unstable. The fifth and sixth eigenvalues then collide at B = 23.39 but 
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FIG. 17: Bifurcation diagram for a stationary anisotropic rod subject to coat hanger boundary 
conditions (oj = 0, R = 0.5512). Bifurcating solutions are not exact helices. 

shortly after become imaginary again in a reverse Hamiltonian-Hopf bifurcation at B — 
23.76. The first and second eigenvalues then also become imaginary again at B = 54.79, 
restabilising the solution (as far as we checked no other eigenvalues cause instability, see 
Fig. [TH1 where A 4 , A 5 and A7 are shown). 

Fig. [20] shows loci of Hamiltonian-Hopf bifurcations in the R-B parameter plane. At 
the end points of the curves the Hamiltonian-Hopf bifurcation coalesces with the pitchfork 
bifurcation. The curve of pitchfork bifurcations is given by ( 193]) and is included in dotted 
lines. We have instability inside the region bounded by solid and dotted curves. Note that 
the stability of the bl solution is almost entirely, but not completely, determined by the first 
two eigenvalues. 

VII. CONCLUSION 

We have shown that whirling current-carrying transversely isotropic rods bifurcate under 
increasing magnetic field (or current) into exact helical shapes provided one applies what 
we call coat hanger boundary conditions. The first bifurcating branch, containing solutions 
with half a helical turn, is stable while higher-order branches are all unstable. 

We stress that the stability analysis of Section 5.1 together with the helical analysis in 
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(a) (b) 




FIG. 18: Evolution with respect to B of the first and second ((a) and (b)) and fifth and sixth ((c) 
and (d)) eigenvalues for an initially stable anisotropic-rod solution on the first bifurcating branch 
(bl). Both pairs of eigenvalues undergo a Hamiltonian-Hopf followed by a reverse Hamiltonian- 
Hopf bifurcation. (7 = 0, oj = 0, R = 0.55.) 

Appendix B gives a complete picture of magnetically-induced helical buckling subject to 
these boundary conditions. In Section 5.1 we show that the critical loads for a straight 
rod are given by equation fl93|) . These bifurcation points are non-degenerate and valid for 
both isotropic and anisotropic rods. The analysis says nothing, however, about the type of 
solutions that bifurcate. They could be helices or not. In Appendix B we then compute 
branches of helical solutions and show that, for isotropic rods (i.e., R — 1), they intersect 
the trivial branch of straight rod solutions precisely at the critical loads computed in Section 
5.1. Taken together these results prove that for an isotropic rod all solutions bifurcating 
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FIG. 19: Evolution with respect to B (7 = and u = 0) of the imaginary part (real part equals 
zero) of the third, fourth and seventh pairs of eigenvalues for non-circular cross-section, R = 0.55. 

from the straight rod are helical. For an anisotropic rod, on the other hand, the critical 
loads are still given by (1931) (with R ^ 1) but no branch of helical solutions in Appendix B is 
found to intersect the trivial branch and therefore the bifurcating solutions are non-helical, 
as confirmed numerically in Section 6.3. 

We have also investigated the stability of post-buckling solutions and found Hopf bifur- 
cations where stable helical solutions lose stability. In the case of a non-rotating anisotropic 
rod we found secondary instabilities given by Hamiltonian-Hopf bifurcations. Unfortunately, 
our method allows us only to study stationary or quasi-stationary (whirling) solutions, so 
it is not clear what type of stable solutions occur after these secondary bifurcations. To 
investigate this one would have to do simulations based on direct discretisation of the PDEs 
([14]) and (US]). 
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FIG. 20: Hamiltonian-Hopf loci for bl solution at u) = 0, R = 0.55. All curves terminate on the 
dotted curve of pitchfork bifurcations given by (|93l) . We have instability of the solution inside 
the regions bounded by solid and dotted lines. Hij labels the curve generated by the ith and jth 
eigenvalue. 

For whirling isotropic rods subject to welded boundary conditions magnetic buckling is 
described by a doubly-degenerate pitchfork bifurcation (see {f}|, where we also showed that 
the same is true for non- rotating states of anisotropic rods). This is because the equations are 
invariant under rotation about the axis {e% axis) of the supports. This symmetry property 
complicated Wolfe's analysis (for non- rotating states of isotropic rods), which had to take 
account of the variational nature of the problem to prove existence of non-trivial bifurcating 
states 5|. By contrast, the coat hanger boundary conditions here introduced break the S 1 
symmetry down to Z2 symmetry (reflection symmetry along Vo and i>i), and no problems in 
the application of standard results from bifurcation theory should arise. Indeed, we find the 
(isolated) critical values of the magnetic field to be given by a remarkably simple explicit 
expression. 

Helical solutions are often used and studied in applications of elastic rods or filaments. 
These solutions are incompatible with the common (aligned) clamped, pinned, Cardan joint 
and other boundary conditions. Boundary conditions are therefore often ignored in these 
studies, making it impossible to do a stability analysis. Here we have introduced a set of 
boundary conditions that does allow for helical solutions, not only in magnetic buckling but 
also in traditional buckling due to compression, twist or whirl. These boundary conditions 
allow one to study exact helical solutions in finite-length rods. We have also proposed a 
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mechanical device that can be used for generating helical solutions in the laboratory. 
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Appendix A: Matrices for the linearisation 



The matrices Bi appearing in equation flbTj) are given by 



( o -k% k\ ^ 



B 1 



.0 







.0 



B, 



( d° u d\ 2 o N 

^21 *^22 
^31 ^32 



F° 



r 3 r 2 

F° 
-F° 



B, 



/ 



-(F 3 )' - 
(F 2 )' - F 3 k° - B(d° 22 dli 



^21^12) 



^21^32) 



(F 3 )' - F?4 

FsKV + Ff^-BidUli-dUiz) 
-(FW-F^ + Bid^-d^) 



-(*?)' 



F?4 



(if?)' - F 2 ^ 



F°k° 



Be 



^ ^11 ^12 ^13 ^ 

^21 ^22 <^23 
y ^31 ^32 ^33 J 



( d° -d° n ^ 



Bp 



2w 



12 
<^22 



\ ^32 ~4l J 

Matrices C, appearing in equation (|68|) are given by 



Ci = B 



C 2 



^ M 3 ° -M 2 ^ 
-M 3 ° Mf 
M 2 ° -Mf i 
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/ /-m /^rl2 r*13 \ 

°3 °3 °3 

^(21 /^22 r<23 

°3 °3 °3 

^31 ^<32 ^33 

°3 °3 °3 J 



where 



C\ l = M°k° 3 + M\k\ - P{uj ■ d° 3 )((f 2 x w • d°) - P(u> • djj)(d£ xw<), 

C 3 12 = (M 3 )' - M°/«2 - PP(<^ • d°)(d° x u ■ d°), 

Cl 3 = -(M 2 )' - Ml4 - Ff + PR(u • d°)(d° x u> • dij) + P(u> ■ d?)(dg x w ■ d°), 

Cf = -(M 3 )' - M 2 °k? + P(w • djj)(d° x « ■ d°), 

Cf = Mini + MX + PR{u> ■ tZ°)(d? xw^) + PP(w • d?)(d° x w • djj), 

Cf = (M°)' - M 2 °/t° - P 2 ° - PP(u; ■ d°)(d? x w • d°) - P(w • d^)(d^ x u - d?), 

Cf = (M 2 )' - M 3 °k? + P° - P(u, • d°)(d° x w ■ d°), 

Cf = -(Ml)' - M$k° 2 + P 2 ° + PR(u ■ d°)(d° x w • d°), 

Cf = M 2 °k° + M°ac$ + P(l - P)(w • d°)(d° x u ■ d°) + P(l - R)(u ■ d°)(d° x uj ■ d§), 



/ 



0-10 
1 




\ 



/ 



1 
OP 
(1 + P) J 



-(d° 2 x w x d?) ■ d? 

P(d° X w x d°) ■ d° 



^ (d§ x w x d°) • d° 
3 = 2P -P(d? x w x d°) ■ d° 

. (d° x w x d°) ■ d° -P(d° x w x d°) ■ d° P(d? x w x d°) • d° - (d° x w x d?) • d 

Matrices -Dj appearing in equation (J7D]1 are given by 
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/ 



-10 
-R 



0- 



r(i+fl) 



Ml 



R 



7 



/ 



V 



-Mi 



T(l+R) 
2 



Rk° 3 
r{i+R) o r(i+fl) o 

2 2 2 1 



-M 2 ° - (1 - R)k' 
Mi - (1 - 



K 2 

-Rk°i 




D A = -jDi. 

All the in the above can be expressed in terms of the moments Mf by means of the 
constitutive relations (G2 
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Appendix B: Helical solutions 



If we assume a helical shape for the magnetically buckled rod then we can derive exact 
solutions as well as buckling loads. For this it is convenient to introduce Euler angles (9, if), 
(p) and relate the director frame {d\, d 2 , d 3 } to the rotating frame {e 1; e 2 , e 3 } as follows: 

d\ = (sin </> cos V> + cos cos 6 1 sin?/') e\ + (sin sin ^ — cos cos 6* cos?/') e 2 
+ cos0sin^e3, 

d 2 = (cos <f> cos if) — sin (f> cos 9 sinif)) e x + (cos <f> sin if> + sine/) cos 9 cos?/') e 2 

— sine/) sin 9 e 3 , (97) 
d 3 = — siia 9 sin if) ei + sin 9 cos ip e 2 + cos 9 e 3 . 

The rate of change of the director frame is given by 

di' = Kxdi (i = 1, 2, 3), (98) 

where k = K±d± + n 2 d 2 + ^3^3 is the curvature vector. Inverting (!98|) and using orthogonality 
of the directors gives 

Ki = d 2 ■ d 3 , k 2 = d 3 ' ■ dx, k 3 = di ■ d 2 , (99) 

which, on inserting fl97|) . yields 

Ki = if)' sin 9 cos (f> — 9' sin <f>, 

k 2 = —if)' sin 9 sin (f> — 9' cos <f>, (100) 
k 3 = <f>' + if' cos 9. 

So for the total curvature k we find 

k 2 = 4 + k 2 2 = 9' 2 + if' 2 sin 2 9. (101) 

Now we make the assumption of a helical centreline and uniformly rotating directors by 
taking 9 to be a constant and setting if)' — Q, <p' — v (Q and v constants). (That <p' may be 
taken constant is not a priori clear. However, performing the following calculation without 
this assumption, while still taking if)' constant, one is quickly led to conclude that <f>' is in 
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fact constant, so to simplify the presentation we assume constancy of 0' from the start.) By 
integrating x' = we then find for the shape 

x(s) = r cos(fis + if)(0)), 

sin^ sin 2 9 

y(s) = r sm(ils + ip{0)), where r = — — — = , (102) 

\l K 

z(s) = L + (s — L) cos 9, 

and with this solution the boundary conditions ( TTHl) - ( |29l) reduce to the four kinematical 
conditions 

cos(fiL + ^(0)) = 0, cos(x + ^(0)) = 0, sin(z/L + 0(0)) = 0, sin 0(0) = 0, (103) 

where we have allowed for a non-zero angle \ between the axes Vq and v\ (see Fig. |3]). 
Compatibility of the director and fixed frame for the case of a straight (but twisted) rod 
(i.e., with 9 = 0) requires 

m+m = l~x, ^L)+ ( p(L) = ^ (104) 

(both frames are aligned at s = L and have a relative rotation x about = k at s = 0). 
Equations ffTUBl and ffTOl) imply 

T17T -\- V T17T 7T 

fi = — ^ = — , 0(O) = O,tt, V(0) = 2 - X - 0(0) (neZ). (105) 

To find the reduced equilibrium equation we integrate the force balance equation (fl6|) to 
get (letting F x = F ■ e 1 , F y = F ■ e 2 , F z = F ■ e 3 ) 

2 * 

F s (s) = -/fi rsin(fis + ^(0)) sin(£2s + ^(0)), 



2 /t 

F„(s) = IB r cos(fis + V(0)) + cos(fis + ^(0)), 



F z (s) = T. 



and insert these expressions, together with (HOOp and the constitutive relations fflQ|) with 
7u = 0, in the moment balance equation (1171) written out in the director frame. The twisting 
moment is directly given by its constitutive relation to be 

M 3 = GJk 3 = GJ(<p' + <0'cos0) = GJ(u + ncos9), (106) 
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i.e., a constant. In the isotropic case (I\ = I2 =: Iq) the moment balance equation for M3 is 
therefore identically satisfied, while the equations for M\ and M2 each give 

(EIo - GJ)n cos 9 - GJu - ^ cos 9 + ^ - cos 9 + P -^- cos 9 = 0, (107) 

Sr it 6 il 

with O and v given by (II 05ft . The torque M of the lower axis Vq about k is given by 

M = M 2 (0) = EI ttsm 2 9 + GJ(v + ttcos9)cos9. (108) 

Thus we have an infinite set of post-buckling helical solution branches parametrised by 
9 G [0,7r/2]. The bifurcating curves in Fig. HJ for instance, are obtained by plotting (nondi- 
mensionalised) r against IB for n = 1,...,6, using equations (11021) and (I107p . The hand- 
edness of helical solutions is determined by the sign of Q, right-handed for Q > 0. The two 
solutions for 0(0) in (11051) correspond to the two branches emanating from the pitchfork 
bifurcations, both with the same handedness but with opposite signs for both K\ and Ki- 

The buckling condition is obtained by setting 9 = 0: 

IB, = ^(n, + X ) 3 - + X f + ^(nn + x ) - B^L + ^ (ri7r + x ), (109) 

or, for x — and in dimensionless parameters, 

uj 2 

B = n 3 7r 3 + nnT + nntu 2 P. (110) 

nil 

This expression shows that an applied tension (T) and inertia (P) stiffen the rod against 
helical buckling, while whirl (uj) softens it, as expected. The critical values of B agree with 
the pitchfork bifurcations in Figs H] and [UJ and with ( 193]) for the present case of an isotropic 
rod. Formula ( I109p is of course also valid in the absence of a magnetic field (B = 0), in 
which case it gives critical values for buckling due to compression (T), twist (x) or whirl 
(co). 

In the anisotropic case (I\ 7^ J 2 ) the moment balance equation for M 3 becomes = 
(EQ 2 + pu 2 )(l2 — I\) sin 2 #sin</>cos(/>, whose only solution compatible with the boundary 
conditions f 1 1 U 5 j) is = or n. But this implies v = and hence, by f 1 1 5 j) . n = and 
Q = x/L- We conclude that in the case x = 0, i.e., in the case of parallel axes v and Vi, 
an anisotropic rod cannot buckle into a helical solution. This agrees with our numerical 
results, which show that the solutions bifurcating at the critical loads ( 193|) when R 7^ 1 are 
non-helical. 
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